x = [ 1 0 0 0 ; 0 1 0 0 ; 0 0 1 0; 1 1 1 0]

x = [x;x]
x = [x;x]
x = [x;x]
x = [x;x]
x = [x;x]

ruido = normrnd(0,1, size(x))

y1 = 10*x + ruido
y2 = 10*x + 2*ruido
y3 = 10*x + 3*ruido

% Componente principal no se modifica
[c1,s1,l1] = princomp(y1);
[c2,s2,l2] = princomp(y2);
[c3,s3,l3] = princomp(y3);

